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TECHNICAL  MILESTONES  PROPOSED  AND  THE  PRESENT  STATUS 

1.  Development  of  Computational  methodology  for  simulation  of  film-condensation:  COMPLETED 

2.  Development  and  testing  of  a  level  set  method  based  code  for  film-condensation  in  2D  Cartesian 
coordinate:  COMPLETED 

3.  Develop  proper  computational  set-up  for  the  present  numerical  simulation:  COMPLETED 

4.  A  detailed  analysis  of  the  results  obtained  by  a  parametric  study  of  the  film-condensation  in  a 
plane  channel:  COMPLETED 

5.  An  attempt  will  be  made  to  develop  another  code  in  3D  cylindrical  coordinate  for  the  film- 
condensation  in  a  pipe:  IN-PROGRESS 


SUMMARY 

The  researchers  in  AFRL  are  trying  to  model  and  analyze  flow  in  a  condenser  of  a  vapor-compression 
refrigeration  system  used  for  thermal  management  of  electronic  components  in  aircrafts.  The  present 
work  is  an  effort  in  that  direction  to  model,  simulate  and  analyze  film  condensation. 

A  novel  dual-grid  level  set  method  based  modeling  and  simulation  of  film-condensation  is 
discussed  here.  An  in-house  code  is  used  for  the  analysis  of  two  different  film  condensation 
problems:  first,  on  the  underside  of  a  horizontal  plate  and  second,  in  a  plane  channel.  The  first  study 
is  done  at  various  degree  of  sub-cooling  of  the  plate,  with  R-123  and  n-pentane  as  the  working  fluid; 
and  various  unsteady  flow  regimes  and  its  effect  on  heat  transfer  is  reported.  The  second  study  is 
done  analytically  for  fully-developed  and  numerically  for  developing  flow  in  a  channel,  with  an 
excellent  agreement  between  the  two  results  in  the  fully  developed  region.  The  simulations  are  done 
for  four  different  cases,  consisting  of  various  standard  thermal  boundary  conditions  on  the  channel 
wall.  The  numerical  results  obtained  for  the  developing  flow  are  consistent  with  the  physical-model 
and  operating-conditions  proposed  in  the  present  work. 

Although  the  present  work  is  fundamental  in  nature,  the  modeling  and  simulation  strategy 
presented  here  can  be  extended  for  to  model  and  analyze  flow  in  a  condenser  for  better  design  and 
integration  in  thermal-management  systems  in  many  aircrafts. 
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TECHNICAL  REPORT 

The  report  consists  of  three  different  sections:  first,  computational-methodology;  second,  code-validation; 
and  third,  computational  fluid  dynamics  study  on  film  condensation.  The  study  is  done  on  two  different 
problems,  film  condensation  on  the  underside  of  a  horizontal  plate  and  in  a  plane  channel. 

1.  Computational  Methodology 

Level-set  method  is  an  Eulerian  technique  to  capture  the  two-phase  interface  (Gada  and  Sharma, 
2009).  The  level-set  function  ( cp )  is  defined  as  a  signed  normal  distance  function  measured  from  the 
interface  such  that  cp  =  0  on  the  interface;  (p>  0  in  the  heavier  fluid,  phase  1  (liquid  phase);  <p<  0  in  the 
lighter  fluid,  phase  2  (vapor  phase). 

Navier-Stokes  equations  are  discretized  by  finite  volume  and  level  set  equations  by  finite 
difference  method.  A  Semi-Explicit  pressure  projection  method  is  used  to  solve  the  NS  equations  on  a 
staggered  grid.Energy  equation  is  solved  using  sharp-interface  method.Temporal  term  is  discretized 
using  Euler  explicit  method  for  NS  equation  and  3rd  order  Range-Kutta  method  for  LS  advection 
equation.  Advection  term  is  discretized  by  2nd  order  TVD  L1NLIN  scheme  for  momentum/energy 
equation  and  5th  order  WENO  scheme  for  LS  advection  equation.  Diffusion  term  is  discretized  by  2nd 
order  central-difference  scheme. 

A  novel  dual-grid  level-set  method  (Gada  and  Sharma,  2011)  is  adopted  in  the  present  study.  It 
utilizes  a  double  refined  grid  to  discretize  the  level-set  function  and  the  temperature,  when  compared 
to  the  pressure  and  velocity  field  variables.  Furthermore,  a  diffused  interface  formulation  is  used  to 
solve  the  flow  field,  while  a  shaip  interface  treatment  is  applied  for  solving  the  energy  equation 
(Gada,  2012).  This  treatment  is  validated  on  a  Stefan  problem,  discussed  in  the  next  section. 

2.  Code-Validation 

Although  the  in-house  code  is  validated  on  Stefan  problem  here,  its  validation  for  film-condensation 
on  the  underside  of  a  horizontal  plate  and  in  a  plane  channel  is  discussed  in  Section  3  below. 

Stefan  problem  represents  a  phase  change  process  progressing  in  one-dimension  over  a  semi- 
infinite  domain  (Alexiades  and  Solomon,  1993),  as  shown  in  Fig.  1.  Initially,  the  entire  domain  is 
filled  with  saturated  vapor.  At  time  r  =  0,  the  left  wall  is  suddenly  sub-cooled  and  maintained  at 
Tx  <Tsat.  This  triggers  a  condensation  process  resulting  in  the  formation  of  a  liquid  layer  between  the 
wall  and  the  vapor.  Thereafter,  the  liquid-vapor  interface  thus  formed  keeps  moving  away  from  the 
wall  due  to  condensation  occurring  at  the  interface,  where  the  vapor  continually  loses  its  latent  heat  to 
the  liquid  side.  At  the  same  time,  a  temperature  gradient  builds  up  on  the  liquid  side  due  to  heat 
diffusion  from  the  interface  to  the  wall  through  the  liquid  layer. 


Sensible  heat  loss  Loss  of  latent  heat 

through  wall  of  condensation 
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Figure  1:  Schematic  of  ID  Stefan  problem 
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This  phenomenon  is  governed  by  an  unsteady  one-dimensional  diffusion  equation  (which  is  solved  in 
the  liquid  phase  only).  In  the  non-dimensional  form,  the  governing  equation  is  given  as 

—  =  -^4-  for  Q<X<XAt)-  r>0  (1) 

dr  dX 2  W 

Here,  Xi  is  the  instantaneous  location  of  the  interface;  6  -  (T  -Tsa,)/(Tsat-  Tw)  is  the  non-dimensional 
temperature,  X  =  x/Lc  is  the  traverse  location,  r  =  ta\ILc2  is  the  non-dimensional  time;  L,  being  the 
characteristic  length  scale,  which  is  taken  equal  to  the  domain  height  for  this  problem. 

The  initial  and  boundary  conditions  are  given  as 

Initial:  XL  (0)  =  0  and  9{X,  0)  =  0  (2) 

Left  boundary :  6*(0,  r)  =  -1  for  r  >  0  (3) 

Right  Boundary:  0{X : ,  r)  =  0  for  z  >  0  (4) 

The  instantaneous  location  of  the  right  boundary,  Xfc),  can  be  found  by  using  Eq.  (5);  which  is 
obtained  by  applying  energy  balance  at  the  interface.  Here,  the  vapor  to  liquid  mass  conversion  rate  at 
the  interface  is  balanced  by  the  heat  diffusing  into  the  liquid  side,  while  neglecting  the  heat  diffusion 
into  the  vapor  side. 

Here,  Sh  =  cpi(Tsat-  Tw)!hn 


dX  j  _  dO 

ij  T  j 

d  r  dX 


for  t  >  0 


(5) 


The  results  obtained  with  n-pentane  as  the  working  fluid  (properties  given  in  Table  1),  on  a  non- 
dimensional  domain  size  of  100x10,  with  a  grid  size  of  100x10  and  Lc  =  0.0001m,  is  compared  with 
the  analytical  results  (Alexiades  and  Solomon,  1993)  and  an  excellent  agreement  is  shown  in  Fig.  2. 


Figure  2:  Comparison  of  interface  evolution  obtained  numerically  and  analytically  for  the  Stefan  problem 


3.  Computational  Fluid  Dynamics  Study  on  Film  Condensation 

For  the  first  time,  condensation  studies  are  done  from  our  in-house  code  (Gada  and  Sharma,  2011), 
developed  for  multiphase  flow.  Two  different  problems  are  simulated  and  analyzed,  discussed  below 
in  separate  subsections. 

3.1.  CFD  Analysis  of  Film  Condensation  on  the  underside  of  a  Horizontal  Plate 

It  is  concluded  from  a  literature  review  that  very  few  studies  have  been  performed  for 
condensation  on  a  downward  facing  plate,  especially  numerically.  The  present  study  aims  at 
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understanding  the  effect  of  increasing  heat  removal  rates  on  the  nature  of  condensate  flow.  With 
increasing  condensate  flow  rates,  the  flow  modes  are  expected  to  undergo  transitions;  as 
suggested  by  studies  on  falling  films  (Dhir  and  Taghavi-Taffeshi,  1981;  Hu  and  Jacobi,  1996). 
The  effect  of  such  a  transition  on  Nusselt  number  is  not  found  in  literature.  Therefore,  the  aim  of 
the  present  study  is  to  elucidate  the  transition  of  condensate  flow  regimes  and  their  effects  with 
increasing  degree  of  sub-cooling. 

3.1.1.  Physical  Description  of  the  Problem 

The  schematic  of  the  domain  for  the  problem  under  consideration  is  shown  in  Fig.  3.  A  flat 
surface,  cooled  from  above,  forms  the  top  boundary  of  the  domain.  Saturated  vapors  lose  heat  to 
this  sub-cooled  surface;  thereby,  forming  a  condensate  film  which  wets  the  underside  of  the  top 
boundary.  The  liquid  film  being  unstable,  under  the  forces  of  gravity  and  surface  tension,  grows 
into  droplets  which  eventually  pinch-off  and  fall  into  the  bulk  vapor.  This  can  be  attributed  to  a 
“Taylor  instability”  that  occurs  at  the  interface  whenever  a  heavier  fluid  sits  on  top  of  a  lighter 
one.  This  instability  grows  in  the  form  of  waves.  The  wavelength  at  which  such  disturbances 
grow  most  rapidly  is  called  the  “most  dangerous  wavelength”  (2d;  given  by  Eq.  (6))  (Bellman 
and  Pennington,  1954;  Sharifi  and  Esmaeeli,  2009).  Because  of  their  rapid  growth  rate, 
disturbances  with  the  most  dangerous  wavelength  are  expected  to  occur  most  often  in  practice. 

Sub-cooled  Wall 

Anti-node  Node  Anti-node 


Outflow  Boundary  Condition 


Figure  3:  Computational  domain  for  condensation  on  the  underside  of  a  horizontal  plate 

3.1.2.  Numerical  Details 

A  domain  of  size  2d  x  2d  is  sufficient  to  capture  the  spatially  periodic  phenomenon  of  droplet-jet 
formation/release  during  condensation.  Similar  domain  definitions  have  proven  to  suffice  in  the 
past  studies  on  film  boiling  (Welch  and  Wilson,  2000;  Tomar  et  al.,  2005)  and  condensation 
(Sharifi  and  Esmaeeli,  2009),  where  the  flow  is  governed  by  the  Taylor  instability.  Furthermore, 
due  to  its  symmetry  along  the  centerline  of  the  domain,  a  computational  domain  ABCD  (shown 
by  the  shaded  region  in  Fig.  3)  of  size  0.52d  x  /.d  suffices  for  the  present  analysis,  with  symmetry 
boundary  conditions  applied  on  its  left  and  right  boundaries. 
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Two  refrigerants  are  selected  for  the  present  study,  namely  R-123  and  n-pentane.  Table  1  shows 
their  properties  and  working  conditions.  A  grid  size  of  100x200  and  125x250  are  found 
sufficient  to  give  grid  independent  results  for  n-pentane  and  R-123,  respectively.  The  liquid  film 
is  perturbed  initially  using  the  curve  given  by  Eq.  (7)  below. 


Ad  =  2 


3a 


g(Pi-Pi) 


(6) 


1 

1 

r  In  x') 

- cos 

40 

64 

v  K d 

(7) 


Table  1:  Properties  of  working  fluids  (  in  SI  units) 


n-pentane 

R-123 

Liquid 

Vapor 

Liquid 

Vapor 

T 

±  sat 

468.3 

455.8 

P sal 

3.3xl06 

3.6x10s 

P 

296.53 

173.22 

686.39 

414.29 

fiX  10s 

3.27 

1.88 

5.52 

3.13 

k 

0.05599 

0.04555 

0.03915 

0.03227 

cpx  1  O'4 

2.39 

3.52 

1.47 

1.91 

a 

3.82xl0‘5 

3.07x1  O'5 

hfg 

5.72xl04 

2.54xl04 

The  boundary  conditions  are  shown  qualitatively  in  Fig.  3.  The  sub-cooled  wall  is  considered  to 
be  cooled  uniformly;  namely  having  a  constant  degree  of  sub-cooling,  A Tsat  =  Tsat  -Tw  =  constant. 
The  boundary  conditions  on  the  left  and  right  boundaries  are  given  by  Eq.  (8)  and  those  on  the 
bottom  by  Eq.  (9).  Top  wall  boundary  conditions  are  summarized  by  Eq.  (10)  below.  The  details 
on  the  equation  and  non-dimensionalization  can  be  found  in  Gada  and  Sharma  (2009).  Elere,  the 
characteristic  temperature  scale  is  taken  as  A Tsat. 


£k=0;  K=o; 

dX  dX  dX  dX 

0;  ^  =  0;P  =  0;  —  =  0;  ^  =  0 
dY  dY  dY  dY 

U  =  0;V  =  0-  —  =  0;  ^-  =  0;  6  =  -\ 
dY  dY 


(8) 

(9) 

(10) 
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3.1.3.  Results  and  Discussions 

Simulations  are  done  for  two  different  working  fluids:  R-123  and  n-pentane  (thermophysical 
properties  and  saturation  condition  given  in  Table  1  above),  at  various  degree  of  sub-cooling  of 
the  horizontal  plate  varying  from  10  to  16  K.  Both  the  working  fluids  indicate  a  transition  in  the 
condensate  flow  regime  from  a  droplet  flow  mode  to  a  jet  flow  mode,  over  a  very  narrow  range 
of  degree  of  sub-cooling. 

Four  flow  regimes  were  clearly  identified  with  increasing  degree  of  wall  sub-cooling,  through 
the  droplet-jet  transition: 

I.  DRIPPING  MODE:  Periodic  single  droplet  release  at  node  and  anti-node  (lOK). 

II.  DRIPPING  MODE:  Periodic  two  droplet  release  at  node  and  one  at  anti-node  ( 12K ). 

III.  DRIP- JET  MODE:  Jetting  mode  at  node  and  periodic  two  droplet  release  at  anti-node  (13K). 

IV.  JETTING  MODE:  Steady  jet  flow  mode  ( 16K ). 

Transition  during  condensation  of  n-pentane 

The  transition  regimes  obtained  for  this  case  are  shown  in  Fig.  4(a)-(d).  Here,  a  series  of 
snapshots  are  shown  to  depict  the  progression  of  flow  over  two  cycles  of  droplet  release,  one  at 
the  node  and  other  at  the  anti-node,  with  increasing  values  of  A Tsat.  Further,  these  can  also  be 
correlated  with  the  instantaneous  surface  averaged  Nusselt  numbers  (Nuavg(z)),  plotted  in  Fig. 
4(e). 

In  regime  I,  Fig.  4(al-a5)  shows  periodically  alternating  single  droplet  release  at  the  node  and 
anti- nodes.  Whereas,  in  regime  II,  Fig.  4(bl-b5)  shows  that  the  number  of  droplets  released  are 
two  at  the  node  and  one  at  the  anti-node.  Therefore,  the  periodic  variation  of  Nuavg( z)  for  the 
release  at  node  and  anti-node  are  same  in  regime  I  (Fig.  4(el))  but  not  in  regime  II  (Fig.  4(e2)). 
Regime  III  experiences  increased  condensate  flow  rates,  which  results  in  the  formation  of  a 
continuous  jet  at  the  node;  while  the  condensate  continues  to  flow  in  the  form  of  drops  at  the 
anti-node  location.  It  is  interesting  to  note  that  each  periodic  droplet  release  cycle  at  the  anti¬ 
node  consists  of  two  droplets.  Therefore,  this  regime  also  attains  a  periodic  variation  of  Nuavg(z); 
with  an  increase  in  its  values  compared  to  the  previous  regimes,  as  seen  in  Fig.  4(e3). 

Finally,  the  flow  transits  to  the  steady  jet  flow  mode  in  regime  IV,  where  the  droplet  release 
at  the  anti- node  continues  initially  in  the  form  of  tip-streaming  (i.e.  drops  are  continuously 
released  with  the  breakup  length  increasing  gradually),  which  finally  leads  to  a  steady  flowing 
jet  formation.  Therefore,  the  corresponding  Nuavg(z )  reaches  a  steady  value.  Such  flow  regimes 
with  increasing  mass  flow  rates,  have  also  been  observed  experimentally  for  falling  films  over 
horizontal  tubes,  as  shown  in  Fig.  5. 
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Figure  4:  Transition  phases  observed  during  condensation  of  n-pentane,  at  various  degree  of  sub 

cooling. 


(a) 

(b) 

■BMBI 

(c)  (d) 


Figure  5:  Flow  modes  obtained  experimentally  during  drop-jet  transition,  with  ethylene  glycol 
flowing  between  horizontal  tube  banks  (a)  droplet  (b)  droplet-jet  (c)  in-line  jet  (d)  staggered  jet 

(Flu  and  Jacobi,  1996) 
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Transition  during  condensation  of  R-123 

The  drop-jet  transition  with  R-123  followed  trends  similar  to  those  obtained  with  n-pentane. 
However,  some  differences  were  seen  in  the  flow  patterns  obtained  in  regime  II  and  III.  In  regime 
II,  a  clear  distinction  from  regime  I  was  not  obtained,  as  the  flow  was  oscillating  between  single 
and  double  droplet  releases  at  the  node  itself.  Nevertheless,  it  must  be  noted  that  the  transition  for 
R-123  happens  over  very  small  temperature  ranges  (see  Fig.  7(b)).  Therefore,  it  is  likely  that 
regime  II  lies  within  a  still  smaller  range;  which  is  not  captured  while  increasing  A Tsat  in  steps. 
Further,  regime  III  was  characterized  by  a  continuous  jet  at  the  node  and  periodic  release  of  a 
single  droplet  at  the  anti-node  (whereas  two  droplet  releases  at  the  anti-node  were  seen  for  n- 
pentane).  This  can  be  attributed  to  the  lower  condensate  flow  rates  obtained  in  the  case  of  R-123. 
This  is  evidenced  from  the  local  plots  of  average  film  velocities  Uavg  and  film  thickness  Afum,  as 
shown  in  Fig.  6.  It  is  seen  that  Um,g.x  is  higher  for  n-pentane  as  compared  to  n-pentane,  while  the 
flow  areas  remain  nearly  identical. 

Surface  averaged  Nusselt  number 

Variation  of  the  time  and  wall  averaged  Nusselt  number  ( Nuavg )  with  increasing  A Tsal  is  plotted 
in  Fig.  7;  and  also  compared  with  a  heat  transfer  correlation  proposed  by  Gerstmann  and  Griffith, 
1967.  Further,  the  progression  of  the  flow  modes  is  also  identified  in  Fig.  7,  with  the  grey  areas 
denoting  erratic  flows  during  transition  between  various  regimes.  Fig.  7  shows  good  agreement 
between  the  present  and  published  results  in  regime  I  and  II  (this  is  expected  since  the  correlation 
has  been  proposed  using  data  obtained  in  the  droplet  flow  mode),  it  disagrees  by  as  much  as  40- 
50%  in  regime  III  and  IV. 


Figure  6:  Local  variation  of  Uavg  and  T/;/,„  obtained  on  the  sub-cooled  wall  with  n-pentane  (R-123)  at 
A Tsat  =  13K  (5K)  and  x  =  318.5  (377) 


Figure  7:  Behavior  of  time  and  wall  averaged  Nusselt  number  for  (a)  n-pentane  and  (b)  R-123 
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Factors  influencing  the  behavior  of  local  and  surface  averaged  Nusselt  number 

In  any  of  the  flow  regimes,  the  heat  transfer  from  the  top  plate  will  depend  on  the  temperature 
gradient  existing  across  the  condensate  film  wetting  the  wall,  the  film  thickness,  the  interfacial 
area  of  the  film  (where  phase  change  occurs)  and  the  average  velocity  in  the  condensate  film.  The 
local  variations  of  some  of  these  variables  over  the  top  wall  are  shown  in  Fig.  8. 

Figure  8  is  drawn  at  time  instants  when  Nuavg(z)  reaches  the  peak  value  in  its  cyclic 
variations  plotted  in  Fig.4(e).  Here,  Nuy  is  the  local  Nusselt  number  evaluated  on  the  sub-cooled 
(top)  wall;  Uavgj  is  the  average  M-velocity  of  the  condensate  film  at  any  location  X;  Ardm.x  is  the 
condensate  film  thickness  at  any  location  X  (therefore,  Af,i,nyX  can  be  considered  to  be  a 
representation  of  the  interface  shape  near  the  wall).  It  is  can  be  seen  that  the  point  of  minimum 
film  thickness  coincides  with  the  point  of  maximum  Nux-  Therefore,  heat  diffusion  through  the 
film  plays  a  dominant  role  in  controlling  the  value  of  Nux- 

Further,  for  both  the  flow  modes,  regime  I  (A Tsat  =  10K)  and  regime  III  (A Tsat  =  13K),  the 
peak  values  of  Nux  are  nearly  identical.  However,  four  peaks  per  wavelength  in  regime  III  are 
responsible  for  a  higher  NuAVG  compared  to  regime  I.  Further,  the  peak  values  of  UavgyX  occur 
slightly  downstream  of  the  point  of  maximum  heat  transfer  coefficient. 


Figure  8:  Local  variation  of  Uavg,  Nu,  Af,im  on  the  sub-cooled  wall  with  n-pentane  (a)  A Tsat  =  10K  at  r  = 

135.0(b)  ATsat=  13K  at  r=  318.5 


3.1.4.  Conclusion 

Direct  numerical  simulations  for  condensation  on  the  underside  of  a  uniformly  cooled  surface 
were  performed  using  a  dual-grid  level-set  methodology.  R-123  and  n-pentane  were  considered 
as  the  working  fluids.  Results  show  that  the  flow  transits  from  a  drop  flow  mode  to  a  steady  jet 
mode  through  the  drop-jet  flow  regimes.  Further,  the  effect  of  the  moving  towards  a  steady  jet 
flow  is  to  enhance  the  average  Nusselt  number  on  the  sub-cooled  surface.  Therefore,  correlations 
available  in  literature  may  under-predict  the  heat  removal  rates  for  the  drop-jet  and  jet  flow 
modes. 
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3.2.  Analytical  and  Numerical  Study  on  Film-Condensation  in  a  Plane  Channel 

In  this  study,  analytical  results  for  fully-developed  and  numerical  results  for  developing  flow 
in  a  plane  channel  with  film-condensation  are  obtained.  However,  the  results  correspond  to 
certain  operating  condition  discussed  below. 

3.2.1.  Physical  Model  of  the  Problem 

Figure  9  shows  the  physical  model,  considered  in  the  present  study,  corresponds  to  hydro- 
dynamically  as  well  as  thermally  developing  vapor-liquid  stratified-flow  in  a  plane  channel;  with 
heating  from  the  top  and  cooling  from  the  bottom  wall.  The  figure  shows  that  the  liquid  and 
vapor  (at  saturated  condition)  enters  from  the  left  boundary  in  fully-developed  stratified  flow 
condition,  with  vapor  flowing  on  top  of  liquid. 

The  flow  configuration  along  with  axial  variation  of  temperature  profile  and  the  two-phase 
interface  is  shown  in  Fig.  9(a),  for  the  channel-walls  subjected  to  different  standard  thermal 
boundary  conditions;  corresponds  to  four  different  cases  (a)-(d)  considered  in  the  present  work. 
Fig  9(b)  shows  the  velocity  profile  along  the  channel  length  for  both  phases.  As  the  interface 
curvature  is  small,  the  surface  tension  force  associated  with  the  interface  is  small  enough  so  that 
the  pressure  jump  across  the  interface  can  be  neglected,  and  both  the  phases  can  be  assumed  to 
remain  at  equilibrium.  As  the  phase  change  takes  place  from  the  interface,  the  interface  is 
assumed  to  remain  at  saturation  temperature  (Tsat). 

The  thermal  BCs  are  applied  in  such  a  manner  that  there  is  always  superheating  from  the  top 
wall  (Tw2>Tsat )  and  sub  cooling  from  the  bottom  wall  (Twi<Tsat).  Up  to  a  certain  length  in  the 
channel  there  exists  a  difference  in  the  wall  heat  flux  between  the  top  and  the  bottom  wall 
(qwiTquU,  which  leads  to  a  jump  in  the  heat  flux  across  the  interface  and  generates  an  interfacial 
mass  flux  from  the  Stephan’s  condition  as  given  in  Gada  and  Sharma  (2009). 

The  region  up  to  which  the  interface  height  changes,  the  phase  change  occurs  along  with 
change  in  thermal  and  velocity  profile;  termed  as  hydro  dynamically  and  thermally  developing 
region.  Beyond  this  region  the  heat  flux  from  the  channel  walls  are  equal  in  magnitude 
(qwl=qw2),  such  that  there  is  no  heat  jump  across  the  interface.  This  leads  to  no  interfacial  mass 
flux  such  that  the  velocity  and  the  temperature  profile  do  not  change  along  the  axial  length.  This 
region  is  called  as  fully  developed  (FD)  region 
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Figure  9:  Physical  model  for  (a)  thermally  and  (b)  hydro-dynamically  developing  flow  in  a 
plane  channel,  with  phase-change  in  the  developing  region.  Developing  (a)  temperature  and  (b) 
velocity  profile  are  shown  in  the  figure.  Arrows  on  channel-wall  in  Fig.  (a)  represents  direction 
of  heat-flux  and  at  the  interface  in  Fig.  (b)  represents  direction  of  interfacial  mass  flow  rate. 
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3.2.2.  Operating  and  Boundary  Conditions 

The  analytical  solution  for  a  condensation  problem  proposed  here  corresponds  to  the  fully 
developed  region,  such  that  the  interface  height  does  not  change  along  the  axial  length. 

Note  that  it  corresponds  to  phase  change  in  the  developing  region  with  no  heat  transfer  in 
the  fully  developed  region,  as  the  heat  gained  from  the  top  becomes  equal  to  that  lost  from 
the  bottom  wall  of  the  channel.  The  thermo-physical  properties  and  boundary  conditions  are 
chosen  such  that  complete  condensation  of  vapor  to  liquid  does  not  occur. 

These  leads  to  certain  operating  conditions,  given  as 


dii!  _ 
~dx  ~ 
dTWl 
dx 


^  =  0,  vi  =  v2  =  o.  ^T  =  0and 

=  dT2  _  dTml  =  dJw2  =  dj2  _  dJm2 

dx  dx  '  dx  dx  dx 


=  0 


(11) 


The  analytical  solutions  are  proposed  here  for  the  four  different  types  of  thermal  boundary 
conditions,  given  in  Table  2. 


Table  2:  Four  different  cases  corresponding  to  various  thermal  boundary  conditions  (CWT: 
constant  wall  temperature;  UHF:  uniform  heat  flux;  and  CHT:  convection  heat  transfer)  on  the 
channel  walls. 


Cases 

(a) 

(b) 

(c) 

(d) 

Top  wall  BC 

CWT:  Tw2>Tsat 

UHF:  qw 

CHT:  h  &  Tx>Tsat 

CHT 

Bottom  Wall  BC 

CWT:  Twl<Tsat 

CWT:  TwI<Tsat 

CWT:  TwI<Tsat 

UHF:  qw 

3.2.3.  Analytical  Solution  from  Energy  Equation 

For  the  operating  conditions  proposed  above,  the  axial  and  transverse  advection  term  in  the 

dT  dT 

energy  equation  are  eliminated  as  =  0  and  v  =  0  respectively.  Thus,  for  each 

phase,  the  steady  state  equation  reduces  to  a  1  -D  steady  state  heat  conduction  equation: 

d2Tf 

k'3^=° . <12> 

where /  =  1  and  2  corresponds  to  liquid  and  vapor  phases  respectively, 

Boundary  conditions  at  the  interface  are  taken  as 
Tf  =  T sat  at y  =  hfd, . (13) 

ki  —  =  k2  —  at  the  interface  . (14) 

1  dy  1  dy  v  ’ 

Integrating  Eq.  (12)  and  using  the  boundary  conditions  at  the  walls  (Table  2)  and  at  the 
interface  Eq.  (13)  &  Eq.  (14),  the  expressions  for  fully-developed  temperature/heat- flux 
profile  (in  the  two-phases)  and  interface-height  are  obtained.  This  is  given  in  Table  3,  for  the 
four  different  cases  considered  in  the  present  work. 
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Table  3:  Analytical  solution  for  fully-developed  two-phase  stratified  flow,  with  condensation  in 
the  developing  region,  for  four  different  cases. 


Cases 

Temperature  Profiles 

Interface  Fleight 

ATi(y) 

AT2(y) 

hfd 

(a) 

( y  Aat 

k±H 

U/d  )  mb'b 

kt  +  k2AT* 

(b) 

Mwt(y-h/d) 

^2 

ki 

Hwt 

(c) 

(  y  \ 

(y  —  byrd)ATsup00 

k1(hH  +  k2) 

J  sub'6 

^2^th,tot 

h{kt  +  fc2A7£) 

(d) 

9wb(y-h/d) 

(y  —  h/d)ATsuPi00 

QwbfaH  +  /^2 )  H-  ^2^^r^sup,oo 

kt 

^2^t/l,tOt 

h-C[wb 

In  the  above  table,  T^y)  =  AT^y)  +  Tsat  ,  T2  (y)  =  AT2  (y)  +  Tsat,  ATsubjb  =  Tsat  - 


ATsup  t  —  Tw2  Tsat ;  ATsup  >co  —  Tsat ,  ATt* 


ATsup  t 
ATsub>b 


AT*  =  ■ 


AT, 


SUp,  00 


AT, 


sub,b 


R 


th.conv 


=  1/h,  R 


th.concL 


—  (.H  —  h fd)/k2,  Rth.tot  —  R 


th.conv 


+  R 


th.concL 


3.2,4.Analytical  Solution  from  Navier  Stokes  Equation 

For  the  fully  developed  conditions  and  the  BCs  at  the  interface,  for  isothermal  flow,  Dutta 
et  al.  (2009)  derived  the  velocity  and  pressure  profiles  from  the  momentum  equation,  This 
is  also  applicable  here  for  the  phase  change,  given  in  Table  4  in  non  dimensional  form. 
Table  4:  Analytical  fully  developed  non-dimensional  velocity  profile  for  each  phases 


Fluid 

UfOO 

1 

A(BY-Y2) 

2 

A(C  +  BY  -  Y2)/q 

A  =  — 


ldP 
2  dX  ’ 


B  = 


dP 

dX 


-12  Qlfd 


1  +  Hfd(rj  —  1) 

1  +  HfdOl  —  1)  ' 

1  +  Hfd  (77  —  1) 


C  = 


.Hk[(?-2Hfd-Hjd)  +  r1Hjd\ 


0?  —  X)Hfd(l  —  HfcL) 

1  +  0?  -  Wfd 

Qlin+xQn 


.  Qlfd 


1  A  XQr,fd 


Qr,fd  — 


(1  -  r))H}d  +  2(317  -  2)H}d  +  3(2  -  3 n)Hfd  +  4(17  -  1  )Hfd  +  1 
17(17  -  1  )Hfd  -  217 Hfd  +  317 Hfd 


n *  _  Ql,in  _  Q 2, 

Vi ,in  —  jj  ’  V2, in  ~ 


UmH  ’ 


uf(Y )  x 

Uf(Y)  =  ^-,X  =  - 


un 


H 


P2  p  =  _PH_  „  =i}fd 
Mi  '  Mit/oo  ’  fd  H 
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3.2.5.Numerical  Solution 

An  in-house  DGLSM  based  code,  (Gada  and  Sharma,  2011)  is  used  for  the  simulation  of 
the  stratified  two  phase  flow  condensation  problem.  The  non-dimensional  height  and  length 
of  the  computational  domain  is  taken  as  1  and  40,  respectively.  At  the  top  and  the  bottom 
wall  of  the  channel,  no  slip  velocity  BC  are  used  and  various  thermal  BC  (corresponding  to 
cases  a  -  d,  given  in  Table  2)  are  used.  At  the  entry  to  the  computational  domain  hydro 
dynamically  fully  developed  boundary  conditions  are  given  (corresponding  to  the  fluids 
viscosity  ratio  rj  and  inlet  flow  rate  ratio  Qr),  using  the  analytical  solution  proposed  by 
Dutta  et  al.(2009).  The  temperature  of  both  the  phases  at  the  entry  are  at  saturated 
condition.  For  the  simulation  Qr  (inlet  flow  rate  ratio)  =  1.272  is  taken. 

The  simulation  is  done  for  steam-water  system  at  near  critical  state,  having  the  state 
parameters  as  psat  =  179  bar  ,  Tsat  =  630  K  and  hfg  =  781  KJ/kg  and  the  thermo  physical 
properties  corresponding  to  the  system  state  are  taken  as 

fluid  1  :  pi  =  538.79  kg/m3 ,  cpi  =  12.6  kJ/kg-K  ,  kj  =  0.41  W/m-K  ,  pi  =  67pN-s/m2 
fluid  2  :  p2  =  133.33  kg/m3 ,  cP2  =  22.1  kJ/kg-K  ,  k2  =  0.13  W/m-K  ,  p2  =  28pN-s/m2 

The  inlet  velocity  to  the  computational  domain  is  taken  to  be  hydro  dynamically  fully 
developed,  whose  interface  height  is  given  by  solving  the  non  linear  equation  in  Hin>  given 
by  Dutta  et  al.  (2009),  as: 

{q(n  -  l)Qr  +  q  -  l)Hin4  +  (4  -  2qQr  -  6q)Hin3  +  (3qQr  +  9q  -  6)Hin2  + 

4q(i  —  q)Hin  —  i  =  o  . (3.1) 


Where  ,  q  =  —  ,  is  the  viscosity  ratio,  Qr  =  , 

M1!  Ql 

Using  Hin  in  expression  of  Uf(Y)  in  Table  4,  the  inlet  velocity  profile  is  obtained.  At  the  exit 
,the  BCs  are  taken  as  fully  developed.  At  the  channel  walls,  the  no  slip  and  no  penetration 
velocity  conditions  are  used.  The  thermal  BCs  for  four  different  cases  are  taken  as  follows: 

Case  (a)  :  Twl  =  620  K,  Tsat  =  630  K,  Tw2  =  655  K 

Case  (b)  :  Twl  =  620  K,  Tsat  =  630  K,  qw2  =  -7  W/m2 

Case  (c)  :  Twl  =  620  K,  Tsat  =  630  K,  T*  =  645  K,  h  =  15  W/m2-K 

Case  (d)  :  qw2  =  -19  W/m2,  Tsat  =  630  K,  T*  =  700  K,  h  =  15  W/m2-K 

The  number  of  grid  points  found  sufficient  for  grid  independent  results  is  800x40,  with 
uniform  grid  spacing  in  X  and  Y  direction  as  AX  =  0.05  and  AY  =  0.025,  respectively. 

3.2.6.Results  and  Discussion 

The  simulation  is  done  for  Rei  =  0.8,  at  near  critical  points  (as  given  before)  in  order  to 
ensure  lower  computational  time. 
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3.2. 6.1.  Fully  Developed  Flow 

Figure  10  shows  an  excellent  agreement  between  the  present  analytical  and  numerical 
results  for  temperature  as  well  as  the  velocity-profile.  The  dashed  lines  in  the  figure 
show  the  analytically  obtained  fully-developed  non-dimensional  interface  height, 
across  which  a  change  in  slope  of  the  profiles  is  seen.  Furthermore,  it  can  be  seen  in 
the  figure  that  the  numerically  obtained  temperature  profile  is  linear  (although  the 
governing  equations  used  for  simulations  included  all  the  terms);  corroborating  our 
proposition  that  the  heat  transfer  in  fully-developed  region  is  only  due  to  transverse 
conduction. 


(d)  U  (c2)  T  (K) 


Figure  10:  Fully  developed  (al-dl)  velocity  (non-dimensional)  profile  and  (a2- 
d2)  temperature  profile  for  all  the  four  cases  (a-d),  along  with  the  analytical 
solutions. 
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3.2. 6.2.  Developing  Flow 

Figure  1 1  shows  the  variation  of  axial  velocity  (non  dimensional)  profile  and 
temperature  profile  for  cases  (a-d).  It  can  be  seen  that  the  result  at  X  =  24  and  X  =  40 
overlaps,  indicating  that  the  flow  becomes  fully  developed  region  beyond  X  =  24. 
Also  the  temperature  plot  shows  non  linear  profiles  near  the  entry  (developing  region) 
and  a  linear  downstream  (developed  region).  This  confirms  the  operating  condition, 
i.e.,  there  is  no  axial  conduction  in  fully  developed  region.  From  the  velocity  plots 
(Fig.  11,  (al)-(dl)),  we  can  see  that  the  volumetric  flow  rate  (i.e,  jUdY,  can  be  seen  as 
area  under  the  U(Y)  curve)  decreases  along  the  axial  direction,  indicating  film- 
condensation  along  the  channel  length. 


>- 


Figure  11:  (al-dl)  axial  velocity  (non-dimensional)  profile  and  (a2-d2)  temperature 
profile  along  the  channel  length,  for  all  the  four  cases  (a-d). 
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Figure  12  shows  the  heat  flux  variation  at  both  the  walls  with  the  axial  length  of  the 
channel  for  all  the  cases  (a-d).  As  the  flow  becomes  fully  developed,  it  is  seen  that  the 
wall  heat  flux  becomes  equal  (qwi  =  qw2),  which  further  justifies  the  linear  nature  of 
the  temperature  profile  as  proposed  in  the  operating  condition. 

Figure  13  shows  the  axial  variation  of  Umax  for  all  the  cases  (a  -  d).  It  is  seen  that 
Umax  decreases  along  the  channel  length,  due  to  condensation  occurring  along  the 
channel  length.  It  is  also  seen  that  the  decrease  in  Umax  is  proportional  to  fully 
developed  interface  height  (which  can  be  correlated  with  1 1  ,cl  given  for  all  cases  in 
other  plots  like  Fig.  10). 


Figure  12:  Axial  variation  of  heat  fluxes  at  the  walls  for  the  cases  (a-d). 


Figure  13:  Axial  variation  of  maximum  velocity  for  the  cases  (a-d). 
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Figure  14  shows  that  the  numerically  obtained  developing  interface  height  increases 
along  the  axial  length  indicating  a  film  condensation  in  the  developing  region.  It 
asymptotes  to  a  fully  developed  value,  which  matches  very  well  with  the  analytical 
results.  Figure  14(d)  shows  a  slight  decrease  in  the  interface  height  near  the  inlet  for 
the  case  (d),  due  to  film-boiling  followed  by  condensation  in  the  downstream  length. 


(c)  X  (d)  X 

Figure  14:  Numerical  obtained  developing  interface,  for  all  the  cases  (a-d).  Fully 
developed  value  compared  with  the  analytical  results. 


3.2.7.  Conclusion 

Analytical  for  fully-developed  and  numerical  for  developing  flow  is  studied  here,  for  film 
condensation  in  a  plane  channel.  Although  the  flow  is  hydro-dynamically  fully  developed  at 
the  inlet,  the  flow  redevelops  due  to  phase  change  in  the  developing  region.  Thus,  analytical 
solution  of  thermally  (thermally  and  hydro-dynamically)  fully-developed  flow  are  obtained. 
The  analytical  solution  is  given  also  for  non-dimensional  interface-height,  volumetric-flow¬ 
rate,  pressure-gradient  and  velocity/shear-stress  profile;  obtained  by  solving  Navier-Stokes 
as  well  as  energy  equations.  Our  analytical  results  can  also  be  used  as  a  test-case,  for 
benchmarking  a  numerical  method  for  two-phase  flow  with  condensation. 

Using  an  in-house  level  set  method  based  code,  2D  simulations  are  done  for  the 
operating  conditions  proposed  in  the  analytical  study.  The  simulations  are  done  for  four 
different  cases,  consisting  of  various  standard  thermal  boundary  conditions  in  the  channel 
wall.  For  the  fully-developed  flow,  an  excellent  agreement  is  found  between  the  numerical 
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and  analytical  results.  The  numerical  results  obtained  for  the  developing  flow  are  consistent 
with  the  physical-model  and  operating-conditions  proposed  in  the  present  work. 

4.  Conclusion 

Out  of  the  five  technical  milestones  proposed  in  this  project  work  (refer  page  no.l),  the  first  four  have 
been  achieved  and  the  detailed  results  are  presented  in  this  report. 
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